String order and adiabatic continuity of Haldane chains and band insulators 
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The ground state of spin-1 Haldane chains is characterized by the so-called string order. We show 
that the same hidden order is also present in ordinary one-dimensional band insulators. We construct 
a family of Hamiltonians which connects adiabatically band insulators to two topologically non- 
equivalent spin models, the Haldane chain and the antiferromagnetic spin-1/2 ladder. We observe 
that the localized spin-1/2 edge-state characteristic of spin-1 chains is smoothly connected to a 
surface-bound state of band insulators and its emergence is not related to any bulk phase transition. 
Furthermore, we show that the string order is absent in any dimensions higher than one. 

PACS numbers: 71.10.Fd, 71.10.Hf, 75.10.Lp, 75.10.Pq 



The antiferromagnetic spin-1 chain is a prominent ex- 
ample of a spin-liquid, an exotic state of matter where 
interactions play a dual role: on the one hand a strong 
Hund's rule coupling locally leads to the formation of 
magnetic spin-1 degrees of freedom, on the other hand 
spin correlations screen magnetism completely. As con- 
jectured by Haldane in the early eighties^, spin-1 chains 
are characterized by a finite gap in their spectrum and 
by exponentially decaying spin-spin correlations. Since 
this seminal work a lot of progress have been made, the- 
oretically and experimentally, towards a full description 
of the so-called "Haldane phase". In particular, Affleck, 
Kennedy, Lieb and Tasaki^ provided a microscopic un- 
derstanding of this correlated ground state, devicing an 
exactly solvable model (AKLT) adiabatically connected 
to the spin-1 chain. The AKLT ground state is made up 
solely of nearest-neighbor singlets (or valence bonds) and 
the gap in the system (and therefore in the spin-1 chain) 
corresponds to the energy needed to break a bond. 

In this paper we want to address the question to what 
extent the Haldane chain is different from an ordinary 
band insulator. The latter has also a gap and an even 
number of electrons per unit cell and, recently, one of 
the authors has shown^ that Luttinger surfaces cannot 
be used to distinguish band- from correlated insulators. 
Two properties, though, are considered peculiar of the 
Haldane phase: First, when a spin-1 chain is cut, a lo- 
calized fractionalized excitation, namely a spin-1/2, is 
formed at the boundary of the chain. These boundary 
spins have for example been observed in the NMR pro- 
file close to the chain ends of Mg-doped Y2BaNiOg^£. 
Second, the ground state possesses a hidden long-range 
order characterizing the entanglement of the spins: Den 
Nijs and Rommelse 6 noticed that, even though true Neel 
order is absent in the ground state, any site with S z = ±1 
is followed by another with S z = =Fl, separated from the 
first by a string of S z = of arbitrary length. This im- 
plies that the so-called string order (SO) parameter 
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A spin-1 chain can be deviced on a spin-1/2 ladder us- 
ing a strong ferromagnetic rung coupling Jr which binds 
two spin 1/2 into a single spin 1. In this case the string 
order parameter takes the formi^i! . 
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and a typical ground state wave function in the Haldane 
phase is shown at the bottom of the first column of Fig. [T] 
As emphasized by Bonesteelii and Kim et. alZ, this wave 
function has a topological feature, namely a vertical line 
(dashed in the figure) crosses always an odd number of 
singlets. As can be easily seen from the figure, this prop- 
erty is directly related to the existence of edge states at 
the boundaries. On the other hand, an antiferromag- 
netic rung coupling leads to another gapful phase, with a 
ground state in a different topological sector of the space 
of singlet wave functions. Here, always an even number 
of singlets is crossed by a vertical and edge states are ab- 
sent (see Fig. [T]). Correspondingl y 7 ! 9 ! 10 the "odd" string 
order S0 o dd vanishes and instead one gets a finite "even" 
string order SO evea 
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is always finite in the Haldane phase. 



In spin-ladder models these two orders never coexist as 
they characterize two different topological sectors. 

In the following we will first show that the ground state 
of an ordinary band insulator can be adiabatically de- 
formed into the ground state of either the Haldane chain 
or the antiferromagnetic ladder. Then we will investigate 
the role both of the string order and of the boundary 
states in the band-insulating phase. Finally, we consider 
the fate of the string order in dimensions larger than one. 

The starting point of our discussion is the following 
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generic Haldane state generic antiferromagnetic state 

FIG. 1: Sketch of the adiabatic path connecting a generic 
Haldane state via a band insulator to a generic antiferromag- 
netic spin ladder. The two ground states in the lower two 
rows belong to the distinct topological sectors of the singlet 
spin Hilbert space as the number of singlet bonds across a 
vertical line (shaded) is either odd or even. 



family of ladder Hamiltonians (see top of Fig. [T]): 

i,a,(T 

+ *« a l,i,(T a 2,i,<r + tDa\ t+l ry a 2 .iM + h.c. 

+ u ^2 n a,iAn a ,%,i + </R^S M S 2 ,i (4) 



where a indicates the two legs of the ladder (or two 
orbitals of an atom). Varying the set of parameters 
[ill t2,tn,tD, U, Jr], it is possible to span a rich phase di- 
agram including a Mott-insulating Haldane phase (U ^> 
hi JR + t R /U < 0), an antiferromagnetic ladder (U ^> t%, 
Jr — 0) an d band-insulators (U = 0, Jr = 0). We 
now describe how these states can be connected adia- 
batically following the route depicted in Fig. [I] First, it 
has been shown both numerically^^ and using analytic 
argument s 7 ' 12 that the Haldane chain can be adiabati- 
cally deformed into a spin-1/2 model with only diagonal 
antiferromagnetic couplings, and, similarly, the antifer- 
romagnetic ladder is connected to a pure rung singlet 
model (second last row of Fig. [I}. These local antifer- 
romagnets can be realized as the large U limit of the 
two-sites Hubbard Hamiltonian 
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for the diagonal or the rung, respectively. For arbitrary 
U /t, the three lowest eigenvalues of Hi oc are a singlet 
with energy E s — —\{U + y/U 2 + 16t 2 ), a triplet with 
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FIG. 2: Energy of the bound state compared to the maxi- 
mum (minimum) of the upper (lower) band as a function of 
tn/tD (ti = —0.1, t 2 — —0.2 and tit = —1). In the inset the 
localization length of the bound state is shown. 



E t = —U and a charge excitation with E c = 
The gap 
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is always finite and, therefore, one can switch off the 
interaction U completely, connecting adiabatically the 
Mott- and the band-insulating ground states. Note that 
the nature of the lowest excitation changes completely 
when E c = E\ while the ground state rotates smoothly 3 
from a state describing two localized spins to one where 
two electrons occupy the bonding band. 

As a final step, we have to show that these two lo- 
cal insulators can be adiabatically connected in the non- 
interacting limit, U, Jr = where the band-structure is 
given by 
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1 ' 2 = (ii + t 2 )cos(fca) ± [(h - i 2 ) 2 cos 2 (fca) + 

(t R - t D ) 2 + 2t R t D (l + cos(fca))] 1 / 2 . (7) 



The states in the second row of Fig. Q] are described by 
ti,ti,tR — and £i,t 2 ,<r> = 0, respectively. For 2\t\_ + 
t 2 \ < \tD + tu\ + \tD ^ tn\ an d h ^ t 2 , the state is 
always a band insulator with a finite gap. Therefore by 
reducing tjj and increasing tu for finite t\ — t 2 we can 
complete the adiabatic path of Fig. Q] For t\ — t 2 the 
Hamiltonian of Eq. (J3]) has a quantum critical point at 
tR = to (related to an extra particle-hole symmetry of 
this model). In Fig. [5] we show that for small t\ — t 2 the 
gap remains finite upon tuning from the diagonal to the 
rung insulator. 

After having shown the adiabatic continuity, we con- 
sider now the evolution of the string order along the adia- 
batic path. In Eqns. |J2]) and ([3]) the S{_ 2 i were operators 
living in the spin Hilbert space while for the Fermionic 
model we have to replace those by Sf_ 2 i = |(ni_2,-i,t — 
ni-2,i,i)- F° r the tw° local singlet states of the third row 
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FIG. 3: Schematic phase diagram of the generalized Hubbard 
ladder Q. The dashed lines denote the adiabatic path of 
Fig. [1] For sufficiently large interactions, either a second or 
first order transition or new intervening phases are expected. 

of Fig. [T] the string order is maximal^, with S'Oodd = \ , 
SO cvcn = and SO cvcn = \, SO od d = 0, respectively. 
Remarkably, the string order is also finite for the two 
local band insulators of the second row: S0 o dd = y^-, 
SO even = and SO cvcn = ^, SO od d = 0, respectively. 
The factor 4 reduction arises because empty and doubly 
occupied states appear with a probability of 1/2, affect- 
ing the spin operators at the beginning and at the end 
of the string. The string, however, still contributes with 
a factor 1. How can these different "order parameters" 
be connected? To answer this question we calculate the 
two string orders along the adiabatic path connecting 
the two local insulators (i.e. on segment B of Fig. [3]). 
With the help of the identities e ac ^ c =: e ^-\)Jc . and 

Sf = -{(: e- 2 <T c M : - : e - 2c ti c ^ :) where : ... : de- 
notes normal ordering, we can reformulate Eq. @ as 

SO odd (i,j) = d: E (e" 2n ^-e- 2 "^) 

r,Z=l,2 
k=j-l 

x exp^ ^2 ( _1 + *) n s,fc,T + (-1 - *)«a,fc,l) 

k=i+l 
s=l,2 

x ( e - 2n '^-T - e- 2ni ^) :) (8) 

and a similar expression is valid for Eq. ^ . In analogy to 
the string order in the spin sector one can also introduce 
a similar quantity in the charge sector, see Appendix [A] 
For a non-interacting system, one can easily eval- 
uate expectation values of the form (: e CiAijCj :) = 
Det[l + G ■ A] using functional integration, where G is 
the equal-time Greens function matrix. We therefore ob- 
tain 

SO* /o (i,i) = i £ (-l) 2 + s ^Det[l+G.A:f a ' k ^], 

a,/3=T4 
r,k=\,1 

(9) 
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FIG. 4: The even and odd string orders, SOcvcn and S'Oodd, 
for the non-interacting ladder as a function of tR/tn (ti = 
—0.1, ti = —0.2, tn — —1). Inset: in the presence of charge 
fluctuations a vertical line can cross both even and odd num- 
ber of singlets. Accordingly the two string orders can coexist. 



where the matrices A e / can be read off from Eq. ([5]) 
and the G matrix can easily be obtained numerically for 
arbitrary ti,t2,tR,tD- Also the determinant in Eq. ([9|) 
can be easily computed for large systems and the SO 
parameters are obtained from the limit \i — j\ 3> £, where 
£ cx 1 /A is the correlation length. 

We find that in one-dimensional band insulators string 
order is always present - even in the absence of inter- 
actions! In contrast to the pure spin models generi- 
cally both string order parameters are finite (see Fig. 3|) 
and they are smooth functions of tn/tp for t\ ^ ti- 
In our model, for tji = to the two string orders 
are equal by symmetry, SO cvcn — S0 o dd (this follows 
from the more general relation SO even (tR,tij,ti,t2) — 
S0 dd{tD,tR,t2,ti) that can be understood in terms of 
the following transformation on the ladder: shift of the 
lower leg by one lattice site and mirror-reflection of the 
ladder respect the x-axis). Only for t\ — ti and tu = to, 
when the gap vanishes, both SO cvcn and S'Oodd deplete to 
because at the critical point singlets of arbitrary length 
are generated. The two orders, S'Oodd and SO cvcn , co- 
exist in gapped insulators also in the presence of strong 
interactions as long as (virtual) charge fluctuations are 
allowed. This can be seen from the inset of Fig. |U a 
doubly occupied or empty site allows to switch from di- 
agonal to rung singlets. 

We now outline the consequences of our discussion 
for the phase diagram of generalized Hubbard ladders, 
sketched in Fig. [3] First, gapped systems are stable 
against small perturbations and therefore there is a fi- 
nite range of parameters which allows to connect adia- 
batically Haldane states and antiferromagnetic ladders. 
Second, in a pure spin model it is not possible to connect 
these two phases as their wave functions live in different 
topological sectors of the spin Hilbert space. Therefore, 
we expect that beyond some critical value of U a phase 
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transition line (first or second order) or an intervening 
phase should separate the even and odd sectors£ i 13 ' 14 i 15 . 

For pure spin models, Starykh and Balents^ 3 - analyzed 
the possible intervening phases and concluded that with- 
out extra fine-tuning there is never a direct second or- 
der transition from the Haldane to the antiferromagnetic 
ladder. Either the transition is first order or intervening 
phases with broken translational invariance appear. The 
absence of a direct second-order transition is consistent 
with our finding that the two high-symmetry phases are 
adiabatically connected. 

To gain further insight in nature of the phase dia- 
gram sketched in Fig. [3] at finite U, we observe that 
for tji = tj) and Jr = 0, our ladder gains an extra 
mirror symmetry (this is evident unfolding the ladder 
on a chain, as shown in the lower part of Fig. O and 
we expect that its properties are similar to the ones of 
the well-studie d 15 ' 16 i 17 ' 18 one-dimensional ionic Hubbard 
model (i.e. a Hubbard model with a staggered poten- 
tial) which has the same symmetries. For this system 
upon increasing U (see upper part of Fig. first the 
charge gap of the band insulator closes and one obtains 
a dimerized phase for U c i < U < U c %. For U > U C 2 
the system becomes a uniform, gappless Mott insulator 
(effectively a uniform Heiscnberg chain 18 ). We therefore 
expect that U c in the schematic phase diagram of Fig. [3] 
has a similar role as U c \, i.e. the charge gap vanishes 
and a dimerized phase (not shown in Fig. [3]) appears for 
U > U c i, tn ~ to, where S'Oovcn is bigger or smaller 
then SO dd, dependingly on which of the two degener- 
ate ground states one considers. According to the results 
of Starykh and Balents 13 discussed above, we do, how- 
ever, not expect a gapless phase for generic parameters 
(e.g. Jr > 0). The role of string order in the charge 
sector of the ionic Hubbard model is briefly discussed in 
Appendix [X] 

At this point we would like to investigate the emer- 
gence of the edge-states following again the path depicted 
in Fig. [T] The localized spin 1/2 at the boundary of the 
Haldane phase evolves smoothly in a surface bound state 
of the "diagonal" band insulator. For finite interactions 
and a proper choice of the chemical potential, this surface 
bound state will be singly occupied giving rise to a local- 
ized spin. In Fig. [2] we show how this edge state merges 
with one of the bands when tn is increased for t\ ^ t^. At 
this point, the localization length of this state diverges 
with {tn/tr> — const.) -1 (inset of Fig. [5]). Nevertheless, 
this is not related to any bulk phase transition. 

Up to now we have only discussed one-dimensional lad- 
ders and it is interesting to study the role of string order 
in higher dimensions. Especially in gapped systems one 
would naively expect that any type of long-range order 
is stable against small perturbations like a coupling to 
neighboring ladders. Using the methods described above, 
we can easily calculate the string order for two- and three- 
dimensional band insulators numerically. Surprisingly we 
find that for arbitrarily weak inter-ladders coupling, t±, 
the string order decays exponentially, SO(i,j) ~ e ~ a \ l ~i\ 
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FIG. 5: Top: the three different phases of the ionic Hubbard 
model as a function of U . A c and A 3 are the charge and 
spin gaps and D = jEj.j" 1 ) 3 ^!,^.' + c] >a Cj+i, a ) is 
the dimerization order parametei^ii^. Bottom: the ladder of 
Eq. unfolded on a chain. An extra mirror symmetry is 
present for = to and Jr — 0. 




FIG. 6: The odd string order of Eq. © for a two-dimensional 
lattice decays as S0 o dd ~ e - "' 1- ^' (with t\ = —0.1, ti = 
-0.2, t R = -1, t D = -2 and t ± = 0, -0.2, -0.4, -0.6 from 
up to down). Inset: a as a function of t\. 



with a ex t\ (see Fig. EJ. A similar observation has been 
made previously by Todo et ah—, who studied spin-1 
ladders numerically. The decay of the string order is a 
consequence of (rare) inter-ladders singlets which intro- 
duce a random phase e ±m l 2 in the string. We note that 
obvious generalizations^ of the string order parameter to 
higher dimensions (replacing strings by squares or cubes) 
are ineffective again due to " dangling singlets" at the in- 
finite surface of such structures. 

In conclusions, we have shown that three phases with 
both spin- and charge gaps and two electrons per unit 
cell, namely Haldane chains, band- insulators and anti- 
ferromagnetic ladders, are actually all the same in the 
sense that their unique ground states can be adiabat- 
ically connected (the excitation spectra, however, may 
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change) . This cousinship is reflected in the prevalence of 
string order in all these one-dimensional phases, even in 
non-interacting band insulators. However, string order is 
not robust in higher dimensions, where it is destroyed by 
arbitrarily weak coupling^. An interesting open question 
is the precise role of the Schricffer- Wolff transformation, 
a central tool to derive effective models. Such a uni- 
tary transformation eliminates charge degrees of freedom 
completely and therefore maps the non-magnetic ground 
state of Hamiltonians with charge gap into one of the 
topological sectors of the spin models. As we have con- 
nected two different sectors adiabatically along our path, 
this seems to imply that the Schrieffer- Wolff transforma- 
tion has to break down at some point. This is consistent 
with a scenario where at U c in Fig. [3] the charge gap 
closes. However, one also has to take into account that 
a Schrieffer- Wolff transformation is not uniquely defined 
as it depends e.g. on the choice of a single-particle basis 
and therefore the topological sector may not be uniquely 
defined in the presence of charge fluctuations. 

The rapid progress in the control of fermionic atoms in 
optical lattices will make possible to study the adiabatic 
evolution of correlated ground states also experimentally 
in the near futur o 21 ' 22 . 

We acknowledge useful discussions with R. Moessner, 
E. Miiller-Hartmann, A. A. Nersesyan, A. Schadschnci- 
der, A.M. Tsvelik and J. Zaanen, J. Zittartz and, espe- 
cially, G.I. Japaridze. We thank for financial support of 
the DFG under SFB 608. 



APPENDIX A: CHARGE STRING ORDER IN 
BAND- AND MOTT-INSULATORS 

In close analogy with the case of the spins (Eqns. @ 
and ([3]) in the text), we can construct string order pa- 
rameters also in the charge sector. We define SO^ vcn and 
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\i—j\^oo \ \ 2 — ' 

l=i+l 

+Sn 2 ,ij (dn ld + 5n 2 ,j)^ (Al) 

i-i 

lim ((($ra M+ i + (5n 2 ,;)exp(i- V" Sn u+1 

|i-j|->oo \ \ A 

l=i+l 

+<Wz) (8nij + i + 5n 2 ,j)} (A2) 
where bni = + n^.j — (rt,;^ + Tfyj,). 
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As shown in Fig. [7j both SO^ dd and SO^ vcn are finite 
for a generic band insulator and - up to a trivial factor 1 /4 
- very similar (but not identical) to the spin string order 
parameters of Fig. 0J In the purely diagonal (t^ = 0) 
and vertical (to = 0) limits only one of the two string 
orders is different from 0. 
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FIG. 7: SO^ dd and SO^. vcn as a function of tii/to for the set 
of parameters t\ = —0.1, ti = —0.2 and tn — —1. 



The presence of charge string order comes as a nat- 
ural consequence of the symmetry existing in the one- 
dimensional band insulator between the spin and the 
charge sectors. In bosonization language, both the 
spin and charge fields are pinned due to the presence 
of a relevant cosine term in the low energy Bosonic 
Hamiltonian^, leading therefore to finite string order 



parameter; 
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As already discussed in the main text, 



this scenario is stable also in presence of interactions (at 
least for truly one-dimensional systems), as long as both 
the Bosonic fields are massive. In the limit of large U, 
the charge string orders will be be strongly suppressed 
by powers of 1/U due to the smallness of the charge fluc- 
tuations (becoming strictly zero for a pure spin model, 
i.e. U — oo) but will remain finite for any finite U. 

In the case of the ionic Hubbard model (corresponding 
to the higher symmetry manifold tn = tp and Jr = in 
the parameter space of Eq. the charge field becomes 
massless at U = U c \ due to the competition between 
two relevant cosine perturbations 15 (the on-site energy 
and the Coulomb interaction). As both for U < U c \ 
and U > U c i the charge field is locked, the charge string 
order will be finite but it will vanish at the gapless point 
U = U c i (the spin string order stays finite at U c \). 
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